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Abstract 

Deep inelastic electron-photon scattering is studied in the range from 1.2 to 30 GeV^ 
using the LEPl data taken with the ALEPH, L3 and OPAL detectors at centre-of-mass 
energies close to the mass of the Z boson. Distributions of the measured hadronic final 
state are corrected to the hadron level and compared to the predictions of the HERWIG 
and PHOJET Monte Carlo models. For large regions in most of the distributions stud- 
ied the results of the different experiments agree with one another. However, significant 
differences are found between the data and the models. Therefore the combined LEP 
data serve as an important input to improve on the Monte Carlo models. 
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1 Introduction 



The measurement of the hadronic structure function F2 crucially depends on the ac- 
curate description of the hadronic final state by Monte Carlo models. The available 
models do not properly account for all features observed in the data, and therefore, 
at present, the accuracy of the measurement of F2 is mainly limited by the imperfect 
description of the hadronic final state by the Monte Carlo models. In previous analyses 
of the individual LEP experiments 1 it had been shown that there are discrepancies 
in several distributions of the hadronic final state between the various QCD models 
and the data. It has also be seen that the data are precise enough to further constrain 
the models. The purpose of this paper is to combine the ALEPH, L3 and OPAL data 
to establish a consistent and significant measurement, which can be used to optimise 
the models. 

In this paper the reaction e+e" — >■ e"'"e~7*7 e+e" hadrons, proceeding via 
the exchange of two photons, is studied in the single tag configuration, where one 
scattered electron^ is detected. The differential cross-section for the deep inelastic 
electron-photon scattering reaction, e(k)7(p) e(k') 7(p) 7*(q) ^ eX, where the 
terms in brackets denote the four-momentum vectors of the particles, is proportional 
to F2{x,Q'^) 0. Here Q'^ = —q^ = —{k — k'Y and x = Q'^/2p ■ q. Experimentally, in 
the single tag configuration, the value of x is obtained using 

... 

where W"^ is the hadronic invariant mass squared, and = —p^ is neglected in calcu- 
lating X. 

The hadronic structure function F2 receives contributions both from the point-like 
part and from the hadron-like part of the photon structure. The point-like part can be 
calculated in perturbative QCD. At low the hadron-like part is usually modelled 
based on the Vector Meson Dominance model. The combined contributions are evolved 
using the DGLAP evolution equation. 

Combining the results of three of the LEP experiments not only reduces the statis- 
tical errors compared to the individual results, but the difference between the results 
also gives a reliable estimate of the systematic, detector dependent, errors. The ex- 
perimental data are fully corrected for trigger inefficiencies and background has been 
subtracted. The experimental distributions presented are also corrected for detector 
effects using different Monte Carlo models, and can directly be compared to the model 
predictions based on generated quantities only, i.e., without the simulation of the de- 
tector response, provided that a well defined set of hadron level cuts defined below 
is applied. This experimental information serves as a basis for improvements on the 
models. 

A set of variables is chosen to compare the corrected data to the hadron level 
predictions of the Monte Carlo models. The variables used are: 

• the reconstructed invariant hadronic mass, PKes; within a restricted range in polar 
angles with respect to the beam axis, 

• the transverse energy out of the plane defined by the beam direction and the direction 

■^Electrons and positrons are referred to as electrons. 
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of the tagged electron, ^'t.out; 

• the number of charged particles, A^trk, 

• the transverse momenta of charged particles with respect to the beam axis, Pt.trk? 

• and the hadronic energy flow, l/A^-d-E/dr], as a function of the pseudorapidity rj = 
— ln(tan(^/2)) with respect to the beam axis, where denotes the number of events. 
The complete definition of how these variables are calculated is given in Section ^. 

For Wj-es, -E-t.out and N^rk the result is a differential event cross-section, which means, 
the distributions have one entry per event, whereas for pt,trk the distribution is a one- 
particle inclusive cross-section. The hadronic energy flow is shown as an average energy 
flow per event, J2E/N, where the sum runs over all objects and over all events in a 
given bin of pseudorapidity. 

The analysis presented here is based on data of the individual experiments taken at 
centre-of-mass energies close to the mass of the Z boson. The approximate ranges in 
and X of the different experiments used in this analysis, and the integrated luminosities 
are listed in Table |I[ The ranges are calculated from the requirements on the energy 
and angle of the deeply inelastically scattered electron, and the x ranges are derived 
from Eq. ^ using the range in and the approximate reach in hadronic invariant 
mass of 3 < < 35 GeV. 





ALEPH 


L3 


OPAL 


range 
X range 


1.2-30 GeV^ 
0.001-0.77 

144 pb-i 


1.2-6.3 GeV^ 
0.001-0.4 

140 pb-i 


1.2-30 GeV^ 
0.001-0.77 

87 pb-i 



Table 1: The approximate ranges in and x of the different experiments and the 
integrated luminosities Cint used in this analysis. 



2 Monte Carlo Models 

The two Monte Carlo models studied are HERWIG5.9 |[ and PHOJET 1.05c 0, using 
the leading order GRV |Q parton distribution functions. 

In the HERWIG model the hard interaction is simulated as eq ^ eq scattering, 
where the incoming quark is generated according to a set of parton distribution func- 
tions of the photon. The incoming quark is subject to an initial state parton shower 
which contains the 7 — q q vertex, and the outgoing partons undergo final state parton 
showers as in the case of e"'"e~ annihilations. The hadronisation is based on the cluster 
model. The initial state parton shower is designed in such a way that the hardest emis- 
sion is matched to the sum of the matrix elements for the resolved processes, g — > qq, 
q — s> qg and the point-like 7 — > qq process. The parton shower uses an angular evolu- 
tion parameter, and so it obeys angular ordering. For point-like events the transverse 
momentum of the partons with respect to the direction of the incoming photon is given 
by perturbation theory. In contrast, for hadron-like events, the photon remnant gets a 
transverse momentum kt with respect to the direction of the incoming photon, where 
originally the transverse momentum was generated from a Gaussian distribution. 
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The HERWIG version used for event simulation is HERWIG5.9+/ct, where the label 
kt denotes that the distribution has been altered from the program default. This 
change is made to improve the agreement for the high-Q^ region between the ALEPH 
and OPAL data and the original HERWIG prediction, where the original program 
is denoted as HERWIG default. The default Gaussian behaviour is replaced by a 
power-law function of the form dk^/ikl + /cq) with k^ = 0.66 GeV. The change is 
motivated by the observation made in photoproduction studies at HERA [1^, where 



the power-law function gave a better description of the data. This is a good example 
of how information from two different, but related reactions can be used to improve on 
a general purpose Monte Carlo program. It is interesting to note that the same value 
of ko is chosen to describe both, photoproduction and deep inelastic electron-photon 
scattering events. The upper limit of k^ in HERWIG-|-fct is fixed at max = 25 GeV^, 
which is almost the upper limit of for the region studied, as shown in Table 0. 

Based on the distributions presented in Section ^ the fixed cut-off has been replaced 
by dynamically adjusting the upper limit on kt on an event by event basis |Tl|. In this 
scheme the maximum transverse momentum k^^^^ is set to the hardest virtuality scale 
in the event, which is of order Q^. The distributions produced with this procedure are 
denoted as HERWIG+fct (dyn). 



The PHOJET Monte Carlo is based on the Dual Parton Model [|T2|. It is designed 
for hadron-hadron, photon-hadron and photon-photon collisions, where originally only 
real or quasi-real photons were considered. It has recently been extended to simulate 
the deep inelastic electron-photon scattering case, where one of the photons is highly 
virtual. For the case of deep inelastic scattering the program is not based on the DIS 
formula using F2 , but the 7*7 cross-section is calculated from the 77 cross-section by 
extrapolating in on the basis of the Generalised Vector Dominance model. The 
events are generated from soft and hard partonic processes, where a cut-off of 2.5 GeV 
on the transverse momentum of the scattered partons in the photon-photon centre- 
of-mass system is used to separate the two classes of events. The hard processes are 
sub-divided into direct processes, where the photon as a whole takes part in the hard 
interactions, and resolved processes. In resolved processes either one or both photons 
fluctuate into a hadronic state, and a quark or gluon of one hadronic state interacts 
either with the other photon, or with a quark or gluon of the second hadronic state. 
Also virtual photons can interact as resolved states, however, the parton distribution 
functions of the photons are suppressed as a function of the photon virtualities. The 
sum of the processes is matched to the deep inelastic scattering cross-section. Initial 
state parton showers are simulated with a backward evolution algorithm using the 
transverse momentum as evolution scale. Final state parton showers are generated with 
the Lund code JETSET [|13| . Both satisfy angular ordering implied by coherence effects. 
The hadronisation is based on the Lund string model as implemented in JETSET. 



3 Experimental Method 

Large data sets are generated with the PHOJET and HERWIG+fct programs respec- 
tively, for ■\/s = Mz- All unstable particles with lifetimes of less than 1 ns are allowed 
to decay in the event generation. In this way the particles of the final state correspond 
approximately to those actually seen in the detectors. The corresponding integrated 
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luminosities for the PHOJET and HERWIG+fct samples are 831 pb~^ and 683 pb~^ 
respectively^. 

The definitions of the phase space and observables include cuts at generator level 
both on the events and on the particles within the events. The cuts are chosen such 
that the individual detectors have good acceptance and therefore detector related un- 
certainties are expected to be small. To select the events at this stage, the following 
cuts are applied to the generated hadron level quantities: 

1. The energy of the deeply inelastically scattered electron has to be larger than 
35 GeV. 

2. The polar angle 6'tag of the deeply inelastically scattered electron with respect 
to either beam direction has to be in the ranges 27 — 55 mrad (low-Q^ region) 
or 60 — 120 mrad (high-Q^ region). The two ranges in scattering angles studied 
correspond to ranges of about 1.2 < < 6.3 GeV^ and 5.7 < < 30 GeV^. 

3. The events are required to contain no electron with energy of more than 35 GeV 
and polar angle above 25 mrad with respect to the beam direction in the hemi- 
sphere opposite to the one containing the deeply inelastically scattered electron. 

4. The number of charged particles A^trk, calculated by summing over all charged 
particles which have a transverse momentum with respect to the beam axis of 
more than 200 MeV and polar angles, ^, with respect to the beam axis, in the 
range 20 < < 160°, has to be greater than or equal to 3. The values chosen 
closely resemble the acceptance of the tracking detectors. 

5. The invariant mass VFres, calculated by summing over all charged and neutral 
particles fulfilling > 200 MeV and 20 < < 160°, corresponding to \'!]\ < 1.735, 
has to be larger than 3 GeV. 

This set of particles and cuts defines the hadron level and the data are corrected 
to this level. The hadron level predictions of W^^s and -Et,out are calculated from 
the charged and neutral particles defined above, and the distributions of A^trk and 
Pt,trk from the charged particles alone. The only exception is the hadronic energy 
flow. For the hadronic energy flow, the same event selection has been applied, but 
all particles are included in the distribution without applying a cut on transverse 
momentum. Figure D shows the differential cross-section dcr/dQ^ within the cuts listed 
above and corrected for detector effects. The vertical line roughly separates the low- 

and high-Q^ regions. Since depends on the energy and angle of the deeply 
inelasticafly scattered electron there is a slight overlap in Q^, but due to the cut on 
^tag the two samples are statistically independent. In the kinematic region studied the 
cross-section prediction of HERWIG-|-fct is about 40% higher than the prediction based 
on PHOJET as shown in Figure 0. 

To study the experimentally observable distributions at the detector level samples of 
60k HERWlG+/ct and 120k PHOJET events, which are statistically independent from 
the samples mentioned above, are passed through the detector simulation programs of 
the ALEPH, L3 and OPAL collaborations, ensuring that all experiments use identical 
events. The objects reconstructed after the detector simulation are energy clusters, 
measured in the electromagnetic and hadronic calorimeters, and tracks, measured in 

''The hadronisation parameters used for the PHOJET and HERWIG+fct simulation have been 
determined from hadronic decays of the Z boson by the L3 and OPAL ||T^ collaborations respec- 
tively. 



4 



the tracking devices. Identical event selection cuts at the detector level are applied 
by all experiments, closely reflecting the cuts applied at the hadron level as described 
above: 

1. A cluster of at least 35 GeV is required in one of the small angle electromagnetic 
luminosity monitors. 

2. The polar angle with respect to either beam direction of the cluster has to be in 
the range from 27 — 55 mrad or 60 — 120 mrad. 

3. The most energetic cluster in the hemisphere opposite to the tagged electron has 
to have energy less than 35 GeV. 

4. At least 3 tracks, fulfilling a set of quality criteria, are required to be observed in 
the tracking devices with 20° < 9 < 160°, and with pt.trk of at least 200 MeV. 

5. The invariant mass, PKes; calculated from all tracks and clusters with pt > 
200 MeV and 20 < 6' < 160° is required to be greater than 3 GeV. 

These objects define the detector level. They are used for the W^es, Et^out, A^trk and 
Pt,trk distributions. The only exception is the energy flow l/A^-d£'/d?7, where again 
no cut on the transverse momentum has been applied. 

With this strategy it is ensured that the hadron level distributions, which are ob- 
tained from the large size samples without detector simulation, and the detector level 
distributions, which are obtained from the samples of smaller size with a detailed de- 
tector simulation for each individual experiment, are statistically independent. Both 
samples are used in the correction procedure applied to the data described in the next 
section. 



4 Corrections for Detector Effects 

Before a measured quantity can be compared to theoretical predictions or to the results 
of other measurements it must first be corrected for various detector related effects, 
such as geometrical acceptance, detector inefficiency and resolution, decays, particle 
interactions with the material of the detector and the effects of the event and track 
selections. Figure^shows the energy flow, l/N-dE/drj, predicted by the HERWIG+fct 
model at the hadron level as well as at the detector level, for the three detectors. The 
events are entered in the figure such that the deeply inelastically scattered electron is 
always at negative rapidities, but not shown. Also shown in the figure is the coverage 
in 7] of the various sub-detectors used in this analysis. A detailed description of the 
ALEPH, L3 and OPAL detectors can be found in [|T^], |T^, and respectively. The 



central region of all the detectors with —1.735 < r] < 1.735 is covered by tracking and 
calorimetry. The forward rj > 1.735, and backward rj < —1.735 regions are covered by 
electromagnetic calorimeters for luminosity measurements and typically extend out to 
7] = ±4.3. The electromagnetic calorimeter of L3 between 1.735 < \ri\ < 3.411 is not 
used in the present analysis. 

For the analysis presented here, the correction of the data to the hadron level is 
done with multiplicative factors, /, relating the measured value X'^^^^ of a quantity X, 
such as a bin content, to the corrected value, X^°^^ using the relation: 

V 

-ycorr -ymeas r -ymeas /q\ 

data data J data vmeas ' V / 
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For distributions, the correction factors are computed bin by bin, e.g. for the energy 
flow, / is the ratio of the hadron level (lightly shaded) and the detector level (darkly 
shaded) distributions shown in Figure ^. In this way of correcting the data the as- 
sumption is made that within the restricted angular range there is little smearing of 
the variables between bins, hence a simple correction factor is justified, and therefore 
no attempt to use an unfolding procedure has been made. The Monte Carlo was used 
to verify the accuracy of this assumption. Application of this correction results in 
measurements corrected to a well-defined kinematical region and particle composition, 
as defined in Section ^. 

The correction factors for the energy flow in the low-Q^ region are shown in Fig- 
ure 1^. The correction factors are near one in the central region of pseudorapidity where 
identical cuts have been applied and they are similar for the three experiments. How- 
ever, there is a much larger spread in the region of larger pseudorapidity r] > 1.735, 
where the experiments have different sub-detectors, different angular coverage, and 
apply different cuts. For ALEPH and OPAL the clusters in the forward detectors are 
required to have an energy of at least 1 GeV, while for L3 this requirement is at least 

4 GeV. This leads to larger correction factors for L3 in that region. In addition, for 
L3, the clusters measured in the forward detectors on the side of the tagged electron, 
i.e. at rj < —2, are not considered in this study. Another difference in the treatment 
of clusters measured in the forward detectors is that OPAL uses a correction function, 
obtained from Monte Carlo studies of the detector response to hadrons, to correct for 
losses in the measurement of the hadronic energy in the forward detector, thereby also 
reducing the correction factors. It should be stressed that for most of the distribu- 
tions used for the comparisons in the following section, only the central detector part 
for polar angles with respect to the beam axis in the range 20 < 6 < 160°, that is 
\ri\ < 1.735, is exploited. The only exception is the energy flow where no cuts on the 
angle are applied. 

The correction factors in the low-Q^ region for the other chosen variables and for 
the ALEPH, L3 and OPAL experiments using the HERWIG+fct and PHOJET mod- 
els, are shown in Figures | and ^ The quoted errors of the correction factors are 
the combined statistical uncertainties of the generated hadron level and the simu- 
lated detector level quantities. While the differences between the correction factors 
obtained from HERWIG-|-fct and PHOJET are very small for the energy flow, they 
can vary significantly for other variables. For example, in the case of OPAL, for W^es 
the HERWIG-|-/ct correction factors are on average about 20% higher than the factors 
obtained with PHOJET. The correction factors for the low-Q^ and high-Q^ regions 
typically differ by less than 10%. 

5 Corrected Data Comparisons 

The discussion of the comparison is subdivided into three parts. First the corrected 
data from the individual experiments are compared to each other and to the Monte 
Carlo models. In this comparison only statistical errors are used and no attempt 
has been made to obtain estimates of systematic errors for the individual experiments. 
Based on the above comparison a modified version of the HERWIG-t-fct model has been 
developed which is described next. Finally the data are combined and compared to the 
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Monte Carlo models. In the combination of the data the spread of the experiments is 
used as an estimate of the systematic uncertainty of the measured distributions. The 
numerical results are listed in Tables p|-[l2| and can be obtained in electronic form [|19|. 



5.1 Data from individual experiments 

The Figures show the corrected differential cross-sections, calculated in the kine- 
matical range described in Section ^, for the data compared with the HERWIG+Zct 
and PHOJET predictions. Figures ^ and |1^ show the data on a linear and log- 
arithmic scale, corrected with the HERWIG+/ct model, while Figures |^, ^ |Tl] and [T3 
show the same data corrected with PHOJET. For all distributions the errors shown are 
the quadratic sum of the statistical errors of the measured quantity and the statistical 
errors of the correction factors. As an example. Table shows the results for the Wj-es 
distribution for the three experiments, listing the statistical error on the data and the 
statistical error on the correction factors, /, separately. 

The experimental results from ALEPH, L3 and OPAL agree with each other within 
errors for large regions in most of the variables studied. However, there are also re- 
gions where they significantly differ from each other, for example, in the region of 
Wres < 10 GeV, Et^out < 5 GeV, for low charged multiplicities and low pt,trk in the 
low-Q^ region. The level of agreement also depends on the Monte Carlo model used 
for correction. The agreement between the experiments is slightly better for the data 
corrected with PHOJET, than for the data corrected with HERWIG-|-fct. In the com- 
bination of the data, the differences between the measured distributions of the different 
experiments will be used as an estimate of the systematic error. 

There are significant differences between the Monte Carlo distributions and the 
data, particularly in the low-Q^ region. The PHOJET distributions lie consistently 
below HERWIG+fct, especially at the low end of the distributions, while the agreement 
of the tails in the high-Q^ region is reasonable. In general the PHOJET predictions 
agree reasonably well with the data, corrected with PHOJET, for large values of the 
variables. However, there are large differences in the low-Q^ region between the data, 
corrected with HERWIG-|-fct, and the HERWIG+/ct predictions in all the distributions. 
For the VTres, -Et,out and Pt,trk distributions (except in the region of low values, where 
the data are inconsistent) the differences between the experiments are much less than 
the HERWIG+fct -PHOJET differences. 

Figures and ^ show the corrected energy flow as a function of pseudorapidity 
in the low-Q^ and high-Q^ region for the individual experiments compared to the 
HERWIG+Zct and PHOJET predictions. Following the energy flow analyses at HERA, 
the statistical errors for the energy flow are taken as a/Z) E'^/N. In Figure [l^the data 
are corrected with the HERWIG-|-/ct niodel and in Figure |l^ with the PHOJET model. 
In the central region of the detector, between —1.5 < rj < 1.5, the data lie between the 
HERWIG-|-/ct and the PHOJET predictions in the low-Q^ region, whereas in the high- 

region the two MC predictions lie closer to each other than the data of OPAL and 
ALEPH. These corrected energy flows in the central region are stable against changes 
in the requirements for the hadronic final state, e.g. changes to the quality cuts of the 
accepted tracks and clusters. 

For rj > 1.5 the data from the various experiments vary much more than the 
statistical errors. Furthermore, it has been found that the corrected energy flows 
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are not very stable against variations of the event selection cuts such as the anti-tag 
condition or a maximum energy cut on the energy deposited in the forward detectors. 
The observed changes are much larger at the detector level than at the hadron level, 
leading to the conclusion that the modelling of the energy response to the hadronic 
energy in the forward region is poor. This can be understood since the sub-detectors 
covering this region have poor energy resolution for hadronic energies and no particle 
identification. Therefore the uncertainty of the hadronic energy flow in the forward 
region is larger than is indicated by the spread of the experiments and it is difficult 
to draw firm conclusions on the description by the Monte Carlo models. However, the 
data appear to lie consistently below the prediction of either generator. 



5.2 Modification of the HERWIG model 

As discussed above, the distribution of transverse momentum of the photon remnant 
with respect to the direction of the incoming photon has been altered, motivated by the 
findings in photoproduction at HERA. At LEP, the modification was initially studied 
as a possible improvement of the agreement between the HERWIG prediction and 
the high-Q2 data of OPAL and ALEPH. While BEKWlG+kt seems to reasonably 
describe the data in the high-Q^ region for low and high values of -E't,out; it overestimates 
the distribution for medium values of -Et,out, as shown in Figure ||. Even though the 
description of the energy fiow is improved with the HERWIG-|-/ct generator, it fails to 
accurately describe the data over a wide range of the pseudorapidity rj. 

While the fixed limit of /c^ ^ax — 25 GeV^ is adequate for the high-Q^ region, in 
the low-Q^ region it introduces too much transverse momentum, which is most clearly 



seen in the transverse momentum distribution of the tracks in Figure |T^. Therefore 
the dynamical adjustment, HERWIG-|-/ct (dyn), as discussed in Section ||, has been 
introduced. As will be seen in the next section, when comparing with the combined 
LEP data, this change leads to an improved description of the data also for the \ow-Q^ 
region. 



5.3 Combined data 

In this section the results from the individual experiments discussed in Section are 
combined. In large ranges of the phase space the individual results agree within the 
statistical precision quoted, however there are also significant differences as discussed 
above. These differences are expected because the previous analyses of the individual 
LEP experiments showed that the systematic errors, which are not included above, 
dominate. Because the Monte Carlo models do not sufficiently well resemble the data, 
evaluating the experimental systematic errors of the measurements based on these 
models will not be very reliable. On the other hand a combined result is desirable 
in order to serve as a guidance for the model builders to improve on their Monte 
Carlo programs. Therefore, in a first attempt to make a combined measurement, 
the experimental systematic error is taken from the spread of the measured results, 
wherever they are significantly different on the basis of the statistical error alone. In 
this case the purely statistical error of the combined result is inflated as discussed 
below. 

The combined distributions from ALEPH, L3 and OPAL are shown in comparison 
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to PHOJET and various predictions from HERWIG in Figures 16-2^. The measured 
values are listed in Tables p|-[T^. The combination of the data follows the procedure 



recommended by the Particle Data Group in section 4.2.2 of [20|. Since this is a crucial 



point of the analysis and because specific choices have to be made in the combination, 
the procedure is briefly discussed below. 

The combined value for a given bin is calculated as the weighted average of the 
measurements of the individual experiments using the statistical errors to calculate 
both the weights and the error of the combined value. To calculate the average bin 
content x and its error cr^ from the individual contents Xi and their statistical errors 
(Ji the following procedure is applied. The average content is calculated from x = 
J2'Wi ■ Xi/J2'^i: using the weights Wi = I /erf- The sum runs over all experiments 
i = 1, . . . , A^ex with non-zero entries in that bin, and the error is taken to be 
The of the average is calculated from = ~ ^i)'^- If in the tail of 

a distribution, for example, some experiments measure zero in a particular bin, then 
X and ^ire scaled by the ratio of the number of experiments with nonzero entries to 
the number of experiments being averaged for that bin. If is larger than A'^ex — 1, 



the error is increased by a factor S = yx^/(-^ex — 1), thereby taking the spread of 
the experiments as an estimate of the experimental systematic uncertainty. Finally, 
to obtain the errors quoted the uncertainties due to the correction factors are added 
in quadrature. Since the same data sets are used by each experiment to calculate 
the correction factors, the corresponding errors are strongly correlated between the 
experiments. To take this into account, the error on the correction factors is included 
by taking the smallest quoted error of the individual experiments as an estimate of this 
systematic error, which is assumed to be 100% correlated amongst the experiments. 

Also in the combination of the results of the individual experiments the l/A^-d-E/d?7 
distribution is treated slightly differently from the other distributions. There is a large 
scatter in the measured values of the different experiments, especially in the forward 
region, as can be seen from Figures and [1^. As a consequence there is also a large 
scatter in the scale factors listed in Tables and [12|. To avoid the combination proce- 
dure manufacturing artificially small errors for bins where the measurements happen to 
coincide, the scale factor applied to obtain the combined measurement is taken as the 
average of the individual scale factors from three neighbouring bins centered around 
the bin under study. 

It is apparent from Figures |I8|, |2g, g2] and g, that the new HERWIG+fct (dyn) 



with the dynamical cut-off lies much closer to the low-Q^ averaged data than the version 
of HERWIG+A;t using the fixed cut-off. However, the difference between HERWIG-|-/ct 
(dyn) and HERWIG default is rather small. 

In the case of the energy flow none of the various HERWIG models is able to 
accurately describe the data. This suggests that even though the new HERWIG+/ct 
(dyn) better describes most of the data distributions, the energy flow is still not well 
understood. 

The PHOJET model. Figures 17, 19, 21, 23 and 25, describes the data reasonably 
well in both regions of Q^, but underestimates the cross-section near the lower limit of 
the distributions. This is understood as a consequence of the high transverse momen- 
tum cut-off of 2.5 GeV for the scattered partons in the hard scattering matrix element. 
Below this cut-off only soft events are generated. Due to the different dependences 
of the soft and hard components in PHOJET, this leads to a suppression of the low-VI^ 
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events 121 



The energy flow corrected with HERWIG+fct and with PHOJET Figures |23 and p5 



mostly agree with each other within errors, except in a few bins in the forward region. 
In this region the data corrected with PHOJET lie below the data corrected with 
HERWIG+fct in the region of the peak at 77 ~ 2 and above the flow corrected with 
HERWIG+fct in the region oi i] > 2.5. 



6 Conclusion 

For the first time the results of deep inelastic electron-photon scattering from three 
of the LEP experiments have been combined and compared to predictions from the 
PHOJET and the HERWIG models. It is found that the data from the ALEPH, L3 
and OPAL experiments agree within statistical errors except near the edges of the 
distributions. Where the spread is larger than expected from the statistical errors, as, 
for example, for low charged multiplicities, this difference is taken as an estimate of 
the detector dependent systematic uncertainty of the measurement. 

In the comparison of the data with the HERWIG+fct model the most striking 
discrepancy is seen in the distributions of the low-Q^ region, where the HERWIG+fct 
model systematically overestimates the data. This discrepancy is found to be mainly 
due to the fixed cut-off for the intrinsic transverse momentum of the quarks in the 
photon in the HERWIG-|-/ct model. By dynamically adjusting the cut-off according to 
the kinematics of the individual event in the HERWIG+/ct (dyn) model the description 
of the data is significantly improved, particularly in the low-Q^ region. 

The PHOJET model describes the data reasonably well in both regions of Q^, but 
underestimate the cross-section near the lower limit of the distributions, due to the high 
transverse momentum cut-off for the scattered partons in the hard scattering matrix 
element. 

The energy flow of the data lies between the predictions of the HERWIG and 
PHOJET Monte Carlo models in the central regions of the detectors. In the forward 
region the Monte Carlo predictions lie systematically above the data. It should be 
noted, however, that it is difficult to assess the systematic errors in this region be- 
cause of the poor resolution of the hadronic energy measured in the electromagnetic 
luminosity monitors. 

The method of combining the data of several of the LEP experiments has proven 
useful to detect shortcomings of Monte Carlo models in the description of these data. 
For the HERWIG Monte Carlo this investigation already demonstrated that the changed 
distribution of transverse momentum kt of the photon remnant with respect to the di- 
rection of the incoming photon, the HERWIG+A;t (dyn) model, gives a better descrip- 
tion of the LEP data. As the data distributions are corrected to the hadron level, they 
can be directly compared to the predictions of the Monte Carlo models, without the 
need of detector simulation, and thus can be used more easily by Monte Carlo model 
builders. 
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Table 2: The individual differential cross-section da/dWj-cs in the low-Q^ region for 
the ALEPH, L3 and OPAL data at ^/s = 91 GeV, calculated in the kinematical range 
defined in the text. The data have been corrected with HERWIG+Zcf The first error 
listed is the statistical error on the data only, the second one is the statistical error 
arising from the correction factors, /. 



14 





low-g2 


high-Q^ 


HERWIG+fct 


MC 




Data 




MC 




Data 




T/f/ 


KjiG V 


(da/dW,,,) [pb/GeV] 


S 


(d(7/dW^res) [pb/GeV] 


S 


Q 
O 


— 41 


10 07 

iy.y 1 


20.25 ± 


i.UD4 


n 
Z 


QQT 
OOl 


D.DOi 


D 


00 1 


1 

± 


U.D4D 


o 
z 


^ AO 


A 

4 


— 


14 10 


13.41 ± 


1 0Q7 


A 


146 


5 304 


4 


i O ( 


1 


370 


1 

X 


301 

OKJ X 




— 6 


10.53 


8 


622 ± 


0.792 


3 


628 


3.961 


3 


510 


± 


0.302 


1 


448 


6 


- 7 


7.732 


5 


669 ± 


0.875 


5 


115 


2.979 


2 


575 


± 


0.219 


1 


006 


7 


-8 


5.718 


3 


832 ± 


0.626 


4 


680 


2.411 


1 


764 


± 


0.175 





419 


8 


-9 


4.106 


2 


485 ± 


0.385 


3 


534 


1.682 


1 


444 


± 


0.165 





204 


9- 


- 10 


3.222 


2 


051 ± 


0.172 


1 


028 


1.243 





854 


± 


0.152 


1 


841 


10 


- 11 


2.323 


1 


292 ± 


0.195 


2 


040 


0.915 





658 


± 


0.140 


1 


813 


11 


- 12 


1.707 


1 


053 ± 


0.127 





957 


0.733 





595 


± 


0.154 


1 


537 


12 


- 13 


1.224 





499 ± 


0.129 


2 


570 


0.459 





318 


± 


0.125 


2 


530 


13 


- 14 


0.926 





433 ± 


0.071 





953 


0.345 





231 


± 


0.065 


1 


322 


14 


- 15 


0.653 





270 ± 


0.069 


1 


366 


0.267 





193 


± 


0.059 





967 


15 


- 16 


0.469 





247 ± 


0.087 


1 


726 


0.180 





136 


± 


0.053 


1 


362 


16 


- 17 


0.345 





129 ± 


0.040 


1 


202 


0.139 





114 


± 


0.045 





455 


17 


- 18 


0.259 





098 ± 


0.032 





874 


0.086 





078 


± 


0.075 


1 


995 


18 


- 19 


0.168 





127 ± 


0.051 


1 


023 


0.064 





047 


± 


0.028 





359 


19 


-20 


0.117 





039 ± 


0.029 


2 


134 


0.061 





039 


± 


0.024 





388 


20 


- 21 


0.087 





062 ± 


0.020 


1 


117 


0.029 





013 


± 


0.010 





561 


21 


- 22 


0.057 





Oil ± 


0.009 


1 


287 


0.020 





009 


± 


0.007 





127 



Table 3: The combined differential cross-section (da/dWrcs) calculated in the k;ine- 
matical range defined in the text for the low-Q^ and high-Q^ regions. The data are 
corrected with the HERWIG+A;t model. The HERWIG+A;t prediction (MC) and the 
scale factor S are shown in addition. 
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Table 4: The combined differential cross-section (da/dWrcs) calculated in the k;ine- 
matical range defined in the text for the low-Q^ and high-Q^ regions. The data are 
corrected with the PHOJET model. The PHOJET prediction (MC) and the scale 
factor S are shown in addition. 
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Table 5: The combined differential cross-section (dcr/dE't^out) calculated in the kine- 
matical range defined in the text for the low-Q^ and high-Q^ regions. The data are 
corrected with the HERWIG+A;t model. The HERWIG+A;t prediction (MC) and the 
scale factor S are shown in addition. 
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Table 6: The combined differential cross-section (dcr/dE't^out) calculated in the kine- 
matical range defined in the text for the low-Q^ and high-Q^ regions. The data are 
corrected with the PHOJET model. The PHOJET prediction (MC) and the scale 
factor S are shown in addition. 
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ld.9o 
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ID. lb 
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Z.oZi) 
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Q 
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2.821 ± 0.414 


Q 1 on 
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i.Uoz 


Q 

y 
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1.540 ± 0.237 


9 074 


0.844 


0.845 ± 0.204 


9 079 


If) 


1.121 


0.979 ± 0.199 


1.839 


0.477 


0.512 ± 0.171 


2.0b9 


11 


0.b42 


0.431 ± 0.097 


1.424 


0.248 


0.15b ± 0.04b 


0.229 


12 


0.34b 


0.195 ± 0.0b5 


1.545 


0.140 


0.101 ± 0.050 


2.04b 


13 


0.180 


0.150 ± 0.032 


1.110 


0.07b 


0.047 ± 0.031 


1.490 


14 


0.083 


0.023 ± 0.013 


0.98b 


0.045 


0.073 ± 0.085 


1.203 


15 


0.058 


0.037 ± 0.01b 


0.554 


0.033 


0.003 ± 0.004 


0.500 



Table 7: The combined differential cross-section (dcr/dA^trk) calculated in the k;ine- 
matical range defined in the text for the low-Q^ and high-Q^ regions. The data are 
corrected with the HERWIG+A;t model. The HERWIG+A;t prediction (MC) and the 
scale factor S are shown in addition. 





low-Q2 


high-Q^ 


PHOJET 


MC 


Data 




MC 


Data 




A^trk 


(da/diVtrk) [pb] 


S 


(d(7/d7Vtrk) [pb] 


S 


3 


11.02 


13.44 ± 0.893 


3 


457 


3 


514 


4.594 ± 0.447 


1.710 


4 


13.7b 


15.73 ± 1.022 


3 


55b 


4 


800 


5.392 ± 0.519 


2.3bb 


5 


10.38 


12.7b ± 0.53b 


1 


890 


3 


878 


4.207 ± 0.343 


1.53b 


b 


7.034 


7.931 ± 0.713 


3 


801 


3 


202 


3.537 ± 0.24b 


0.48b 


7 


3.891 


4.405 ± 0.443 


3 


180 


2 


025 


2.031 ± 0.292 


2.2b3 


8 


2.249 


2.412 ± 0.321 


3 


147 


1 


319 


1.385 ± 0.203 


1.881 


9 


1.099 


1.28b ± 0.172 


2 


151 





735 


0.805 ± 0.173 


1.917 


10 


0.558 


0.b84 ± 0.13b 


2 


243 





340 


0.387 ± 0.07b 


0.999 


11 


0.279 


0.318 ± 0.054 


1 


171 





179 


0.158 ± 0.049 


1.2b3 


12 


0.13b 


0.127 ± 0.035 


1 


373 





091 


0.113 ± 0.037 


0.97b 


13 


0.072 


0.053 ± 0.042 


2 


829 





05b 


0.017 ± 0.008 


l.ObO 


14 


0.03b 


0.030 ± 0.017 





990 





024 


0.018 ± 0.011 


0.479 


15 


0.018 


0.029 ± 0.014 


1 


054 





013 


0.003 ± 0.003 


0.500 



Table 8: The combined differential cross-section (da/dA^trk) calculated in the kine- 
matical range defined in the text for the low-Q^ and high-Q^ regions. The data are 
corrected with the PHOJET model. The PHOJET prediction (MC) and the scale 
factor S are shown in addition. 
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iow-g2 


high-g2 


HERWIG+A;t 


MC 


Data 




MC 


Data 




Pt,trk 


(da/dpt,trk) [pb/GeV] 


S 


(da/dpt,trk) [pb/GeV] 


S 


0.2 


-0.3 


757.2 


728.1 


± 


42.83 


6.778 


265.5 


224.7 


± 


17.29 


3.641 


0.3 


-0.4 


656.1 


600.9 


± 


61.32 


10.96 


230.7 


201.3 


± 


6.290 
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0.4 


-0.5 


509.7 


473.4 


± 


36.16 
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± 
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2.412 


0.5 


-0.6 
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± 


26.63 
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146.1 


125.0 


± 


13.71 
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0.6 


-0.7 


292.4 


259.5 


± 


19.57 


5.247 
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± 


4.512 


0.312 


0.7 
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223.5 
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± 


17.54 
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± 
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0.8 
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± 
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69.09 
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± 
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1.411 


0.9 
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± 
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± 
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1.0 
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± 


5.713 


2.696 
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± 
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0.575 
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53.99 


± 
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36.74 
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± 
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± 


2.616 


1.151 


30.71 
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± 


4.872 


3.061 


1.3 
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53.41 


27.85 


± 


2.930 


2.345 


24.91 


20.61 


± 


2.116 


0.283 


1.4 
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44.89 


22.39 


± 


4.662 


4.713 


21.39 


21.69 


± 


2.352 


0.503 
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- 1.6 


36.04 


18.35 


± 


1.716 


1.466 


17.05 
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± 


1.441 
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30.36 


13.68 


± 


1.789 
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14.64 
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± 


1.488 
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24.71 
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± 


1.309 


1.595 


13.07 
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± 


1.403 
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20.61 


9.879 


± 


1.162 
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10.92 
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± 
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16.67 


7.824 


± 


1.246 


1.704 


9.516 
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± 


1.385 
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13.29 


5.296 


± 


0.712 
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7.548 


6.173 


± 
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0.916 
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9.261 


3.743 


± 


0.395 


0.826 


5.498 


5.039 


± 


0.973 
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6.303 


2.671 


± 


0.351 


0.513 


3.843 


2.861 


± 


0.624 


1.740 


2.6 


-2.8 


4.399 


1.865 


± 


0.342 


1.310 


3.082 


2.250 


± 


0.396 


0.096 


2.8 


-3.0 


3.155 


1.248 


± 


0.348 


1.947 


2.365 


1.273 


± 


0.253 


0.306 


3.0 


-3.4 


1.962 


0.778 


± 


0.212 


2.669 


1.471 


1.620 


± 


0.163 


0.982 


3.4 


-3.8 


1.087 


0.617 


± 


0.152 


2.685 


0.868 


0.832 


± 


0.129 


0.785 


3.8 


-4.2 


0.651 


0.516 


± 


0.077 


0.718 


0.524 


0.295 


± 


0.078 


1.555 


4.2 


-4.6 


0.344 


0.234 


± 


0.210 


9.881 


0.245 


0.199 


± 


0.090 


0.000 


4.6 


-5.0 


0.146 


0.101 


± 


0.032 


1.584 


0.128 


0.107 


± 


0.072 


0.000 



Table 9: The combined differential cross-section (dcr/d]?t,trk) calculated in the kine- 
matical range defined in the text for the low-Q^ and high-Q^ regions. The data are 
corrected with the HERWIG+A;t model. The HERWIG+A;t prediction (MC) and the 
scale factor S are shown in addition. 
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± 
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± 
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0.4 
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± 
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± 


9.387 


2.093 
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17.09 
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± 
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15.41 
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± 
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0.839 
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12.58 
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± 
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± 
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2.703 


2.664 


± 
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0.211 
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0.321 


0.359 


2.8 


-3.0 


0.831 


0.752 
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0.137 
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1.536 
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0.602 
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0.545 


0.580 


± 


0.179 


2.801 


0.990 


1.239 


± 


0.428 


3.554 


3.4 


-3.8 


0.337 


0.280 


± 


0.095 


2.388 


0.545 


0.779 


± 


0.141 


1.407 


3.8 


-4.2 


0.211 


0.189 


± 


0.036 


2.066 


0.274 


0.649 


± 


0.093 


0.000 


4.2 


-4.6 


0.132 


0.123 


± 


0.075 


2.943 


0.187 


0.097 


± 


0.235 


0.000 


4.6 


-5.0 


0.084 


0.079 


± 


0.050 


3.143 


0.090 


0.009 


± 


0.065 


0.000 



Table 10: The combined differential cross-section (dcr/dpt,trk) calculated in the kine- 
matical range defined in the text for the low-Q^ and high-Q^ regions. The data are 
corrected with the PHOJET model. The PHOJET prediction (MC) and the scale 
factor S are shown in addition. 
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low-g2 


high-g2 


HERWIG+A;t 


MC 


Data 




MC 


Data 




V 
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dE/dr])[GcV] 


S 


{1/N- 
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-3.0 - 


-2.5 


0.264 


0.281 ± 


0.054 


0.000 
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0.076 ± 


0.066 


2.947 


-2.5 - 


-2.0 


0.605 


1.010 ± 


0.054 


2.364 


0.432 


0.651 =b 


0.108 


2.759 


-2.0- 


-1.5 


1.477 


1.715 ± 


0.078 


0.743 


1.332 


1.449 ± 


0.096 


1.119 


-1.5 - 


-1.0 


2.290 


2.225 ± 


0.056 


1.451 


2.516 


2.221 ± 


0.088 


0.352 


-1.0 - 


-0.5 


2.164 


2.077 ± 


0.043 


0.446 


2.616 


2.366 d= 


0.076 


2.040 


-0.5 - 


-0.0 


1.879 


1.794 ± 


0.033 


0.885 


2.362 


2.365 ± 


0.100 


0.915 


0.0 - 


0.5 


1.891 


1.707 ± 


0.032 


0.554 


2.304 


2.383 ± 


0.135 


3.555 


0.5 - 


1.0 


2.252 


1.993 ± 


0.045 


1.492 


2.533 


2.602 ± 


0.202 


4.606 


1.0- 


1.5 


3.028 


2.559 ± 


0.120 


1.770 


3.169 


3.054 ± 


0.372 


3.627 


1.5 - 


2.0 


3.379 


2.754 ± 


0.226 


4.235 


3.337 


2.879 ± 


0.369 


8.051 


2.0 - 


2.5 


3.100 


1.780 ± 


0.256 


6.335 


3.080 


1.888 ± 


0.321 


1.702 


2.5 - 


3.0 


2.808 


1.063 ± 


0.319 


3.643 


2.735 


1.059 ± 


0.150 


1.018 


3.0- 


3.5 


1.915 


1.267 ± 


0.254 


8.168 


1.918 


1.496 ± 


0.209 


0.825 


3.5 - 


4.0 


1.122 


0.992 ± 


0.249 


4.124 


1.089 


0.528 ± 


0.133 


3.557 


4.0- 


4.5 


0.579 


0.634 ± 


0.144 


2.831 


0.548 


0.292 ± 


0.100 


1.992 



Table 11: The combined energy flow {1/N ■ dE/drj) calculated in the kinematical 
range deflned in the text for the \ow-Q^ and high-Q^ regions. The data are corrected 
with the HERWIG+fct model. The HERWIG+A;t prediction (MC) and the scale factor 
S are shown in addition. 



22 





\ow-Q^ 


high-Q2 


PHOJET 


MC 


Data 




MC 


Data 




V 


(1/iV- 


dE/dr])[GeV] 


S 


(1/iV- 


dE/dr])[GcV] 


S 


-3.0 - 


-2.5 


0.511 


0.248 ± 


0.030 


0.245 


0.176 


0.055 ± 


0.043 


2.748 


-2.5 - 


-2.0 


1.064 


1.218 ± 


0.046 


0.903 


0.578 


0.606 ± 


0.098 


4.332 


-2.0- 


-1.5 


1.836 


1.879 ± 


0.058 


1.028 


1.739 


1.507 ± 


0.117 


0.041 


-1.5 - 


-1.0 


2.170 


2.165 ± 


0.046 


1.457 


2.683 


2.353 ± 


0.109 


2.151 


-1.0 - 


-0.5 


1.835 


1.974 ± 


0.043 


1.576 


2.580 


2.456 ± 


0.107 


3.213 


-0.5 - 


-0.0 


1.601 


1.771 ± 


0.033 


1.268 


2.312 


2.419 ± 


0.131 


1.085 


0.0 - 


0.5 


1.604 


1.730 ± 


0.047 


1.715 


2.260 


2.348 ± 


0.156 


4.633 


0.5 - 


1.0 


1.941 


1.983 ± 


0.063 


3.785 


2.529 


2.429 ± 


0.234 


5.302 


1.0- 


1.5 


2.503 


2.340 ± 


0.139 


2.296 


3.080 


2.523 ± 


0.348 


5.272 


1.5 - 


2.0 


2.770 


2.108 ± 


0.275 


7.846 


3.277 


2.349 ± 


0.371 


7.974 


2.0 - 


2.5 


2.718 


1.707 ± 


0.284 


10.74 


3.109 


1.707 ± 


0.312 


2.817 


2.5 - 


3.0 


2.772 


1.046 ± 


0.337 


2.135 


3.061 


0.934 ± 


0.175 


1.052 


3.0- 


3.5 


2.565 


1.536 ± 


0.187 


6.334 


2.762 


1.446 ± 


0.196 


1.249 


3.5 - 


4.0 


1.929 


1.467 ± 


0.164 


1.052 


2.150 


0.821 ± 


0.207 


2.550 


4.0- 


4.5 


1.193 


1.065 ± 


0.090 


0.393 


1.217 


0.411 ± 


0.223 


3.190 



Table 12: The combined energy flow {1/N ■ dE/drj) calculated in the kinematical 
range defined in the text for the low-Q^ and high-Q^ regions. The data are corrected 
with the PHOJET model. The PHOJET prediction (MC) and the scale factor S are 
shown in addition. 
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Figure 1: The differential cross-section, da/dQ^, observed in the two ranges in scatter- 
ing angles ^tag studied, compared to the predictions of the HERWIG-|-A;t and PHOJET 
models. The cross-sections are given in the kinematical ranges described in the text. 
The errors are statistical only. 



24 




OPAL detector coverage 




2 4 6 

Ti=-log(tan(e/2)) 



ALEPH TPC time projection chamber 
ECAL el. mag, calorimeter 
HCAL hadronic calorimeter 
LCAL luminosity calorimeter 
SICAL silicon tungsten calorimeter 

L3 TEC time expansion chamber 

EE endcop el. mag. calorimeter 
EB barrel el. mag. calorimeter 
ALR lead— scintillator calorimeter 
FECAL forward el. mag. calorimeter 

OPAL CJ central jet chamber 

CAL el. mag. / hadronic calorimeter 
FD forward lead scintillator det. 
SW silicon tungsten calorimeter 



Figure 2: The HERWIG+/ct energy flow, ■ dE/dr], at the hadron level (lightly 
shaded) as well as on detector level (darkly shaded), as measured by the three detectors. 
The coverage of the different detector components are indicated by dot-dashed lines. 
In the case of L3 the ALR is not used in this analysis. 
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PHOJET correction factors 
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2 4 6 

Ti=-log(tan(e/2)) 



Figure 3: The HERWIG+A;t and PHOJET correction factors, /, for the ALEPH, L3 
and OPAL energy flow for the low-Q^ region. The symbols are shghtly displaced for 
better visibility. The vertical lines indicate the central rapidity region, \r]\ < 1.735. 
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Figure 4: The HERWIG+A;t correction factors, /, for tlie ALEPH, L3 and OPAL 
Wres, -St.out, -^trk and Pt,trk distributioHS for the low-Q^ region. The symbols are shghtly 
displaced for better visibility. 
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PHOJET correction factors 
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Figure 5: The PHOJET correction factors, /, for tiie ALEPH, L3 and OPAL H^.es, 
-St.out, -^trk and pt,trk distributions for the low-Q^ region. The symbols are shghtly 
displaced for better visibility. 



28 




Figure 6: The W^es distributions from ALEPH, L3 and OPAL for the low-Q^ (left) 
and high-Q^ region (right), corrected with the HERWIG+Zct niodel on a hnear scale 
(top) and on a log scale (bottom) . 
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Figure 7: The W^^ distributions from ALEPH, L3 and OPAL for the low-Q^ (left) 
and high-Q^ region (right) , corrected with the PHOJET model on a hnear scale (top) 
and on a log scale (bottom) . 
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HERWK; corrected low Q" data 
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Figure 8: The Et,out distributions from ALEPH, L3 and OPAL for the low-Q^ (left) 
and high-Q^ region (right), corrected with the HERWIG+Zct niodel on a hnear scale 
(top) and on a log scale (bottom) . 
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Figure 9: The £;t,out distributions from ALEPH, L3 and OPAL for the low-Q^ (left) 
and high-Q^ region (right), corrected with the PHOJET model on a hnear scale (top) 
and on a log scale (bottom) . 
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HERWK; corrected low Q" data HERWIC corrected high Q" data 




Figure 10: The A^trk distributions from ALEPH, L3 and OPAL for the low-Q^ (left) 
and high-Q^ region (right), corrected with the HERWIG+Zct niodel on a hnear scale 
(top) and on a log scale (bottom) . 
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Figure 11: The A^trk distributions from ALEPH, L3 and OPAL for the low-Q^ (left) 
and high-Q^ region (right), corrected with the PHOJET model on a hnear scale (top) 
and on a log scale (bottom) . 
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HERWK; corrected low Q" data 



^000 
% 900 
> 800 



HERWK; corrected high Q" data 



HERWIG+k, 
PHOJET 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



HERWIG+k, ALEPH 
PHOJET ^ OPAL 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Figure 12: The pt,trk distributions from ALEPH, L3 and OPAL for the low-Q^ (left) 
and high-Q^ region (right), corrected with the HERWIG+Zct niodel on a hnear scale 
(top) and on a log scale (bottom) . 
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Figure 13: The pt,trk distributions from ALEPH, L3 and OPAL for the low-Q^ (left) 
and high-Q^ region (right) , corrected with the PHOJET model on a hnear scale (top) 
and on a log scale (bottom) . 
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HERWIG corrected low Q data 
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Figure 14: The hadronic energy flow from ALEPH, L3 and OPAL for tlie low-Q^ 
(left) and high-Q^ region (right), corrected with the HERWIG+fct model. The shaded 
band indicates the forward region of the experiments. 
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Figure 15: The hadronic energy flow from ALEPH, L3 and OPAL for the low-Q^ 
(left) and high-Q^ region (right), corrected with the PHOJET model. The shaded 
band indicates the forward region of the experiments. 
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HERWIG corrected low Q data 
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Figure 16: The combined Wj-es distribution from ALEPH, L3 and OPAL for the \ow-Q^ 
(left) and high-Q^ region (right), corrected with the HERWIG+fct model on a linear 
scale (top) and on a log scale (bottom). The data are compared to three different 
model assumptions of the HERWIG+fct model. 
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Figure 17: The combined VFres distribution from ALEPH, L3 and OPAL for the low- 
(left) and high-Q^ region (right), corrected with and compared to the PHOJET 
model. 
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HERWIG corrected low Q data 
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Figure 18: The combined -Et.out distribution from ALEPH, L3 and OPAL for the 
low-Q^ (left) and high-Q^ region (right), corrected with the HERWIG+fct model on a 
linear scale (top) and on a log scale (bottom) . The data are compared to three different 
model assumptions of the HERWIG+Zct model. 



PHOJET corrected low Q data 



n 



if 



-|- combined data 
PHOJET 




10 12 14 




^ 10 



s 8 

6 

' 4 
2 




PHOJET corrected high Q data 



-f- combined data 
PHOJET 




10 



12 



14 




10 



12 



14 



Figure 19: The combined -Etjout distribution from ALEPH, L3 and OPAL for the 
low-Q^ (left) and high-Q^ region(right), corrected with and compared to the PHOJET 
model. 
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HERWIG corrected low Q data 
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Figure 20: The combined Ntrk distribution from ALEPH, L3 and OPAL for the low-Q^ 
(left) and high-Q^ region (right), corrected with the HERWIG+fct model on a linear 
scale (top) and on a log scale (bottom). The data are compared to three different 
model assumptions of the HERWIG+Zct model. 
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Figure 21: The combined A^trk distribution from ALEPH, L3 and OPAL for the low-Q^ 
(left) and high-Q^ region (right), corrected with and compared to the PHOJET model. 
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Figure 22: The combined pt,trk distribution from ALEPH, L3 and OPAL for the 
low-Q^ (left) and high-Q^ region (right), corrected with the HERWIG+/ct model on a 
linear scale (top) and on a log scale (bottom) . The data are compared to three different 
model assumptions of the HERWIG+Zct model. 
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Figure 23: The combined pt,trk distribution from ALEPH, L3 and OPAL for the low- 
(left) and high-Q^ region (right), corrected with and compared to the PHOJET 
model. 
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Figure 24: The combined hadronic energy flow from ALEPH, L3 and OPAL for tfie 
low-Q^ (left) and fiigli-Q^ region (right), corrected with the HERWIG+fct model. The 
data are compared to three different model assumptions of the HERWIG+fct model. 
The shaded band indicates the forward region of the experiments. 
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Figure 25: The combined hadronic energy flow from ALEPH, L3 and OPAL for the 
low-Q^ (left) and high-Q^ region (right), corrected with the PHOJET model. The data 
are compared to the PHOJET model. The shaded band indicates the forward region 
of the experiments. 
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